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The accompanying Comment by A. Holas and N. H. March 
[Phys. Rev. A 66, 066501 (2002)] is concerned with the issue 
of whether or not kinetic energy can be represented by an ef- 
fective local potential, as required for an exact Thomas- Fermi 
theory equivalent to Kohn-Sham density-functional theory. 
They dispute [R.K. Nesbet, Phys. Rev. A 65, 010502(R) 
(2001)], which concludes that for more than two electrons 
the use by Kohn and Sham of the Schrodinger kinetic energy 
operator is variationally correct, while the equivalent local 
potential required for a valid Thomas-Fermi theory, a Frechet 
functional derivative of the Kohn-Sham ground-state kinetic 
energy functional, does not exist. The argument of Holas and 
March is clearly invalid for the simple example of the low- 
est triplet state of a two-electron atom with noninteracting 
electrons. Why this fails, as do earlier arguments in the liter- 
ature, has been explained in recent publications, summarized 
here. 

The principal issue in dispute is the question whether 
an exact ab initio Thomas-Fermi theory (TFT) exists 
for ground states, equivalent to Kohn-Sham density- 
functional theory (DFT). This depends on the definition 
of a ground-state kinetic energy functional of the elec- 
tron density and on the definition and computability of 
its density-functional derivative. Holas and March claim 
that a local (Frechet) functional derivative of the Kohn- 
Sham ground-state kinetic energy functional exists, as 
a rigorous consequence of variational theory. This can- 
not be true, because it would imply the equivalence of 
ground-state TFT and DFT, which result in clearly in- 
consistent variational equations for systems with more 
than two electrons [1]. 

The conclusions of the cited paper [2] are misstated 
by Holas and March [3]. In contrast to their allusion to 
"... Nesbet's worrisome claims...", which they hold to be 
contrary to the accepted opinion that "...DFT is funda- 
mentally correct and internally consistent...", the anal- 
ysis summarized in [2] concludes that Kohn-Sham the- 
ory correctly uses the kinetic energy operator t = — ^V 2 
of Schrodinger, and is variationally correct in the local- 
density approximation(LDA). However, there is no for- 
mal proof that a local (Frechet) functional derivative ex- 
ists in general for the exchange-correlation density func- 
tional, nor for the Kohn-Sham kinetic energy functional. 
In fact, the latter Frechet derivative cannot exist in a 
valid quantum theory for more than two electrons. This 
precludes an exact TFT for ground states. Because Holas 
and March cite standard arguments which lead to simi- 
lar inconsistencies, it is important to trace the point of 



failure of these arguments. This will be done here. 

The fundamental incompatibility of Thomas-Fermi 
theory (TFT) and Kohn-Sham theory (DFT) for more 
than two electrons has not been widely recognized. If the 
argument of Holas and March and earlier conclusions in 
accepted literature were sound, TFT would be correct. If 
the two theories are incompatible, then DFT must fail, 
in sharp contrast to the history of applications of these 
theories. To illustrate this incompatibility, consider the 
lowest 3 S state of an atom with two noninteracting elec- 
trons. Kohn-Sham theory solves the pair of Schrodinger 
equations 

i(f>2s = {t2s - v(r)}<t)2 S . (1) 

The orbital wave functions are orthogonal because the 
energy eigenvalues are unequal. They are independently 
normalized to unity and are equally occupied. An equiv- 
alent TFT would use a single Lagrange multiplier p, for 
the global normalization constraint J pd 3 r = 2. If a local 
density functional derivative of the kinetic energy func- 
tional T exists in the form /Spj = Vt(y), this implies the 
Thomas-Fermi equation 

Vt(t) = fJ, - v(r), (2) 

Eqs.(l) cannot be deduced from Eq.(2), because two in- 
dependent Lagrange multipliers cannot be derived from 
one. Unconstrained by the exclusion principle, this 
Thomas-Fermi equation has a ground-state solution in 
which the partial density p\ s is normalized to 2, while pis 
is normalized to zero. The total energy, 2ei s , is below the 
correct ground state. The theory fails because a single 
global Lagrange multiplier p does not suffice to normalize 
pis and pis separately. It cannot replace the set of inde- 
pendent Lagrange multipliers c t required for Fermi-Dirac 
electrons. Any theory that produces a Frechet derivative 
of a kinetic energy functional implies incorrect physics 
unless all e; are equal. 

Holas and March present an argument which, if cor- 
rect, would support the conclusion that solution of a 
Thomas- Fermi equation and of Kohn-Sham equations for 
the same variational model must lead to the same results. 
This is obviously false for the simple model cited above. 
Eqs.(9) of Holas and March give the functional differ- 
ential ST of the defining Schrodinger orbital functional 
T = Because p = J2i n iPi = in 

DFT, a sum of orbital densities, 

6T = J2m f d 3 r{a -v(r)}8 Pl (r) 
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for variations about solutions of the nonintcracting 
Schrodinger equations. This is valid for all density varia- 
tions induced by orbital variations, whether or not the 
orbitals are normalized [2]. These equations demon- 
strate that partial variations SiT of the orbital func- 
tional T about stationary states are unique functionals 
of the partial density variations 5pi [4] . This implies that 
the Kohn-Sham ground-state functional T s [p] can be ex- 
tended to include such infinitesimal variations and defines 
a density functional derivative for each orbital subshcll 
[2,4], 



This is a Gateaux functional derivative, the generaliza- 
tion to functional analysis of a partial derivative [5]. A 
Frechet derivative, equivalent to a local potential func- 
tion vt, exists only if all are equal [4]. 

The Kohn-Sham density functional T s is obtained by 
evaluating T for ground-state orbital functions, restrict- 
ing functional T s to normalized ground states of other- 
wise arbitrary external local potential functions. The 
fallacy in the argument of Holas and March is that they 
cannot deduce the functional derivative from this func- 
tional, limited to normalized orbitals and densities. They 
modify Eqs.(9) by replacing the parameters ej, implied 
by the Schrodinger equation, by the single parameter p, 
appropriate to TFT. This leads to their Eq.(ll), which 
is not consistent with Eqs.(9) unless all integrals J Spi 
vanish, as they do if all pi are separately normalized. 
Since this requires independent Lagrange multipliers 
replacing them by the single parameter p is not justified 
unless all €i are equal. For fixed orbital normalization, all 
integrals of constants drop out, so that Eqs.(9) and (11) 
cannot be distinguished. No conclusion can be drawn 
about functional derivatives unless the definition of T s 
is extended to density variations driven by orbital varia- 
tions that are unconstrained in the orbital Hilbert space. 

Lagrange multiplier formalism requires the relevant 
functional to be defined for variations unconstrained by 
normalization [7] . Otherwise Euler-Lagrange variational 
equations cannot be derived. This motivates the exten- 
sion of T s as defined by Kohn and Sham to infinitesi- 
mal function-neighborhoods of ground-state densities p, 
generated by infinitesimal variations of the orbital func- 
tions <j>i [4]. For nonintcracting electrons, Fermi-Dirac 
statistics requires separate normalization of each elec- 
tronic subshcll density, which implies independent La- 
grange multipliers a indexed by orbital subshclls. Un- 
less all these orbital energies are equal, they cannot be 
replaced by a single parameter p. TFT uses such a global 
parameter to normalize the total density, with no refer- 
ence to shell structure. This is the reason for the historic 
failure of TFT to describe electronic shell structure [2] . 

For stationary states, the Gateaux derivative implied 
by Eqs.(9) is operationally equivalent to the linear oper- 
ator t, because (ej — v)(j>i = t<j)i [2]. It can be expressed 
as a local mean orbital kinetic energy, 



ST _ #(r)f&(r) 
niSpi(r) pi(r) 

For independently normalized partial densities pi, this 
determines modified Thomas-Fermi equations equivalent 
to the usual Kohn-Sham or Schrodinger equations for 
nonintcracting electrons [6]. This restores Fermi-Dirac 
statistics and orbital shell structure to Thomas-Fermi 
theory 

Holas and March cite the work of Englisch and En- 
glisch [8,9] as supporting their argument. This is incor- 
rect, because the rigorous discussion of [8,9] cannot dis- 
tinguish between Frechet and Gateaux functional deriva- 
tives. Only normalized density functions are considered, 
and constants in the Euler-Lagrange equations arc not 
determined. Englisch and Englisch [8] prove the exis- 
tence in general of a Gateaux derivative. They do not 
consider the subshell structure of the DFT density func- 
tion, but arbitrarily insert a non-indexed single parame- 
ter A into their Eq.(4.1), which if correct would imply an 
exact TFT. As shown above for the simple example of the 
lowest ls2s 3 S state of an atom with two noninteracting 
electrons, this result is physically incorrect. Construction 
of a Gateaux derivative consistent with the operator t in 
nonintcracting Schrodinger equations and in Kohn-Sham 
equations uses only the mathematics of variational the- 
ory in the orbital Hilbert space, following a long tradition 
of textbook derivations [7] . 
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